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We use numerical integration to solve the field equations of conformal gravity, assuming a metric 
that is static and spherically symmetric. Our solution is an extension of that found by Mannheim 
and Kazanas [l|; it indicates, as expected, that gravitation in this model should be attractive on 
small scales and repulsive on large ones. 
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I. INTRODUCTION 

Conformal gravity (CG), along with other alternate 
theories of gravity, has been described in a review by 
Mannheim )2|; this paper will referred to as PM. There 
have been important developments since this review was 
written [§}, but they will not concern us here. 

CG is based on the conjecture that at the deepest level 
the laws of nature should be conformally invariant. Grav¬ 
itation theory is therefore based on an action principle 
derived from the Weyl tensor, (PM, sections 5.1 

and 8.7). An early achievement in the theory of CG was 
the discovery by Mannheim and Kazanas of the analytic 
solution of the field equations 

= 0 (1) 

for a static, spherically symmetric source, analogous to 
the Schwarzschild solution in conventional gravitation 
(Mannheim and Kazanas [2j; this paper will be referred 
to as MK. See also PM, section 9.2). The line element 
has the form 

rlr 2 

ds 2 = —B(r)dt 2 H-— + r 2 dUL 2 (2) 

B(r) 

in the notation of Weinberg Q ■ 

MK showed that outside a static, spherically symmet¬ 
ric source the function B(r) is given by 

B(r) = 1 — — ^2L1 _ 3 ^-y -|- yr — kr 2 (3) 

r 

The final kr 2 term is important only at cosmological dis¬ 
tances, and serves to embed the metric in a de Sitter 
universe. The term proportional to 1/r is analogous the 
the Schwarzschild solution. The yr term is new, and has 
been used by Mannheim and collaborators [sj to explain 
galactic rotation curves without recourse to dark matter. 

Mannheim has also constructed a conformally invari¬ 
ant cosmological model (PM, sections 3.5, 8.7, 10). In ad¬ 
dition to the usual fermion fields, he introduces a scalar 
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field, S, that undergoes a symmetry breaking transition 
and acquires a constant non-zero vacuum expectation 
value, Sq. The resulting field equations are (PM, equa¬ 
tion (188)): 


4a g W^ (4) 

where a g is a dimensionless coupling constant, and the 
energy-momentum tensor, T ^ v , is given in PM, equation 
( 66 ): " 


T» v = iip^{x)[d v +T v {x)\it- ^g^hSoW 

~ ^ 0 2 (V - \g^R a a ) (5) 

The fermion contribution will not concern us here, be¬ 
cause we will deal only with the vacuum equations. But 
we note that even if the fermion fields are zero, is 
not zero, but depends on S 0 and R 

T I1V = -^ 0 2 (b>‘ v - \g^R a }j ( 6 ) 

In developing this cosmological model, Mannheim 
points out that for the FRW space normally assumed, 
is identically zero, so the field equations are simply 

= 0 , (7) 

identical to the equations of conventional cosmology but 
with an effective gravitational constant that is negative 
(PM, equation (224)): 


G eff = 


3c 3 

AttSI 


( 8 ) 


II. STRUCTURE OF THE FIELD EQUATIONS 
FOR A STATIC, SPHERICALLY SYMMETRIC 
SOURCE 


We start with the observation that the solution ob¬ 
tained by MK for a static, spherically symmetric source 
did not use the field equation (H]) but the simpler Jl]). 
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Parameter 

Value 

Units 

T surface 

1.0E20 

m 

T max 

1.0E26 

m 

center 

1.0E23 

m 

V 

-1.0E-46 

m -2 

P 

5.39E13 

m 

7 

1.97E-28 

m -1 

k 

(-1.42E-51) 

m -2 


TABLE I: Parameters used in the integration. The second 
column gives values appropriate for the galactic cluster NGC 
3198. v surface is the disc radius, and r max is chosen to be a 
typical cosmological distance, reenter is the geometric mean 
of rsurface and r max . g is set equal to -1/reenter- The value 
for k is enclosed in parentheses because it is not specified at 
the beginning, but is determined during the fitting process. 



FIG. 1: B — 1 vs r, using outward integration with 
log(r/r ce nter) as the independent variable. This graph uses 
parameters for NGC3198, as listed in Table [T] 


Combining with the complete field equations 
can be written 

W ,lv + 77 (V" - ^ v R a a J = 0 (9) 

where 77 = Sq/ (24a s ) is a constant of dimension 
length - . We shall assume the magnitude and sign of 
?7 are completely unknown, and for the purposes of this 
paper can be chosen as we wish. We note, however, that 
some numerical work on the matching of interior and ex¬ 
terior solutions ( 6 ] suggests that 77 is positive. We will 
call the first term in © the “main term”, and the rest 
the “eta terms”. Since 77 is not dimensionless, we can 
expect that the main term will dominate for small r, and 
the eta terms for large r. 

If the eta terms alone are set to zero we get the equa¬ 
tion 

R" v - l -</ LV R° a = 0 (10) 

which leads to the same line element as before, ©, with a 
metric function of the form B(r) = 1+a/r + br 2 , without 
the linear term of MK. 

This leads to the conjecture that at small scales the so¬ 
lution of MK is appropriate, but that at some distance of 
order 77 -1 / 2 the solution goes over to one of Schwarzschild 
form. The main purpose of this paper is to demonstrate, 
by numerical integration, that a solution of that kind 
does, indeed, exist. 


III. THE INTEGRATION: ASSUMPTIONS AND 
PARAMETERS 

As shown in table HI we use parameters approximately 
equal to those given by Mannheim for NGC3198 [h]. Inte¬ 
gration is carried out over six orders of magnitude, from 
t surface to a cosmological distance, r max . 


The line element has the form © at the beginning, 
when the main terms of the field equation dominate, and 
at the end, when the eta terms dominate. In the tran¬ 
sition region we cannot be so sure, and a more exact 
treatment might use the more general line element with 
two unknown functions: 


ds 2 = — B{r)dt 2 + C(r)dr 2 + r 2 dVt 2 (11) 


For simplicity, however, we will assume that, for weak 
fields, the form © is adequate throughout. 

The tensor equation © gives three different equations, 
but because of the relations between them it is only nec¬ 
essary to integrate one. The rr component yields the 
extension of MK (14); primes denote differentiation with 
respect to r: 


B 


-1 


W rr 

1 


77 ( R rr - -g rr R a a 
1 


B'B'" - — ( B"Y - — (BB"’ - B’B") 
6 12 v ’ 3r v ’ 


+ 


1 

3r 2 . 
V 


BB" 
B" + - 


(- b'Y 


2 B 2 1 

3 ?’ 3 3 r 4 3 r 4 


2 B 
r 2 


( 12 ) 


Setting the right side equal to zero we get a third order 
equation that can be integrated straightforwardly. 

The linear term in the solution of MK is sufficiently 
small in practice that B(r ) never deviates far from unity. 
We therefore integrate the equation for the function 
A(r) = B(r) — 1: 


A'" = j-l(A ") 2 + J-(A'A") 


1 

3 r 2 


(A + i) A" + (A'y + —— {A + 1) A' 
J 3r 6 


A{A + 2) 
3 ?’ 4 



(A +1 A'\ 

V 3 r 6 ) 


( 13 ) 
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Some details of the integration are given in the ap¬ 
pendix. Here we simply present figure [1] In the fitting 
process, the parameter k has been varied to make the 
curve reach — 1.0E —6 at r max . We have arranged for this 
small r 1 2 3 4 5 6 7 8 component at the end to match a de Sitter space; 
at the right of the graph, the curve is headed for a hori¬ 
zon at r = 1.0E29, where A = — 1 and B = 0. We note 
that |fc| <C 1771 , a necessary condition if the cosmological 
term is only to become noticeable for r > r cen ter- 

We see, first, that for suitable choices of 77 and fc, a 
solution of the field equations does exist that duplicates 
the solution of MK for small r, and gives a Schwarzschild- 
like solution for large r. Further, we find that 77 should 
be chosen negative. If we change the sign of 77 , we find 
that even for values of \k\ as large as l/r 2 enter the final 
value of A is always large and negative. 

IV. COMMENTS 

Our graph indicates an attractive force for r < r cen ten 
and a repulsive force for r > r cen ter, in accordance with 
the well-known feature of Mannheim’s model, that grav¬ 
itation is apparently attractive on small scales but repul¬ 
sive on large ones. 

Checking figure Q] against observations may prove diffi¬ 
cult, since so much of the figure corresponds to distances 
greater than 1.0E22m. We emphasize that in this paper, 
for simplicity, we have assumed space to be free of matter 
for r > ?" sur face- Mannheim and O’Brien js] have pointed 


out, however, that distant matter can have an observable 
effect on galactic rotation curves, and our analysis would 
have to be extended to take this into account. 

As a first estimate of how distant matter might affect 
galacic rotation curves if figure [l] is appropriate, let us 
consider a simplified model in Minkowski space. Sup¬ 
pose sources within a distance of r cen ter produce a linear 
gravitational potential, kr , with k > 0. Summing over 
all these sources, and approximating the sum by an inte¬ 
gral, we can show straightforwardly that they produce a 
potential near the origin that goes like cir 2 , wth C\ > 0 . 
In other words, a test particle released a short distance 
from the origin will be drawn back towards it. 

Suppose in addition that sources more distant than 
r center generate a potential more suitable for the right half 
of figured] namely a Newtonian potential, M/r, added 
to a cosmological component, Kr , with k < 0. Imagine 
that these sources are distributed uniformly out to some 
horizon, r max . The integrated effect of these sources near 
the origin can be shown to be a potential C 2 r 2 , with C 2 < 
0 . 

We see that local and distant sources have opposite 
effects, and determining which will prevail is a delicate 
matter, requiring a full relativistic treatment, This will 
not be attempted here. We observe, however, that if, af¬ 
ter the S — i Sq transition, there is a time interval during 
which the horizon distance is smaller than r cen ter, then 
all sources will be attractive, and we may regain the type 
of cosmology familiar from Einsteinian relativity. 
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with respect to z, the equation becomes 


A = -2A + 3A + 


4(A + 1) -2A 


-(-A+i) +4 A(-A+A 
(A + 1) (-A + A 


A 2 


+ 8 (A + l)A-4A(A + 2) 

+ 3 ? 7 ro exp(2z) A + A — 2A^j | (Al) 


Appendix A: Details of the integration 


We use the integrator RADAU5, as described in [j] 

For initial conditions we use ©, the solution of MK, 
evaluated at r sur f ace . 

Rather than use equation 01 directly, we transform it 
in two ways. First, we change the independent variable to 
z = log(r/r cen ter)- Using dots to denote differentiations 


The second transformation is simply a scale change. 
A(r ) is always much less than unity, and its derivatives 
are smaller still. We can get more manageable magni¬ 
tudes by multiplying all variables by a constant, P, cho¬ 
sen so that PA(r sul - f ace ) = 1. The scale factor is removed 
after the integration is complete. Writing Ap = PA , 
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Ap = PA, etc., 

Ap = — 2Ap + 3Ap + 


4(A P + P) -2A P _ 

^^—Ap + Ap^ + 4 Ap ^—Ap + Ap^ 


(Ap + P ) Ap + Ap^ + A 

+ 8(Ap + P)A P - 4 Ap(A P + 2P) 

+ 3Pr)rQ exp(2z) Ap + Ap — 2Ap^j | 


(A2) 


(Because A(r) is always so small, we can use either this 


complete equation or its linearized form; the resulting 
graphs are very similar.) 


Before beginning the main integration, we find two val¬ 
ues of k, one of which drives A(r) to a large positive value 
at r max , while the other drives A(r) to a large negative 
value. These two values are used as input to a root¬ 
finding program, such as rtbis from Numerical Recipes 
[§]. Each pass through the root finder involves a complete 
integration, with the output being the value of A(r max ). 
k is automatically adjusted in this way to bring A(r max ) 
to the chosen final value. 







